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The present work reports an experimental observation of thermal entanglement in a clusterized 
spin chain formed in the compound Na2Cu5Si40i4. The presence of entanglement was investigated 
through two measured quantities, an Entanglement Witness and the Entanglement of Formation, 
both derived from the magnetic susceptibility. It was found that pairwise entanglement exists 
below ~ 200 K. Tripartite entanglement was also observed below ~ 240 K. A theoretical study of 
' entanglement evolution as a function of applied field and temperature is also presented. 
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' Since the early years of quantum mechanics, entanglement has attracted much attention due to its fascinating 
^ features, such as non— locality, as exemplified in the Einstein, Podolsky and Rosen paradox Recently, it has been 
discovered that entangled states constitute a valuable resource for quantum information processing and it has 
■ raised a great number of studies about entanglement in many different quantum systems. Over the past few years, 
much effort has been done in developing methods to detect and quantify entanglement. 

Until a few years ago, entanglement was not believed to exist beyond atomic scale. The most common arguments 
against entanglement on larger scales is that macroscopic objects contain a large number of constituents that interact 
^2 ■ with its surroundings, inducing the decoherence phenomena which leads to loss of entanglement as size, complexity 
and system's temperature increases. Surprisingly, it was theoretically demonstrated 0, [j| that entangled states can 
exist in solids at finite temperature and this kind of entanglement is referred in literature as "thermal ent ang lement 
" . Since the publication of several theoretical works jLii^. 1. 1, J> S Ell, El, El, El, El El, E3, El, El, Eol , a few 
experimental evidences have been reported [2l|, [H, IH, [2J| , confirming the presence of entanglement in solids state 
' systems. 

\^ \ The study of entanglement in solid state physics (for a detailed review see [11]) is of great relevance to the area 
\^ • of Quantum Information and Quantum Computation, since many proposals of quantum chips are solid state-based. 
CO I Furthermore, the demonstration that entanglement can change the thermodynamical properties of solids, such as 
magnetic susceptibility [2l|, shows that entanglement can be related to significant macroscopic effects. Hence, this 
T-H ■ subject establishes a interesting connection between quantum information theory and condensed matter physics. 
C , The task of entanglement quantification is still an open problem in general case (for recent reviews see [26l. [27l. [28j ) . 
' Many quantities have been proposed to quantify entanglement, one of them is the Entanglement of Formation (EF) 
[H, HI], which can be easily calculated in the case of two spin 1/2 particles, but it can not be measured directly 
in most of the cases. Hence, usually the detection of entanglement is done using a quantity called Entanglement 
Witness (EW) [3l[, which is an observable that, by definition, has a positive expectation value for separable states 
^ . and negative for some entangled states. 

In this work, we report an experimental observation of thermal entanglement in the spin chain system formed in 
the compound Na2Cu5Si40i4 [H, 111], by using a thermodynamical EW based on the magnetic susceptibility [34| . 
Furthermore, we have also derived, using the complete knowledge of the Hamiltonian of the spin chain, a relationship 
between the Entanglement of Formation of each pair of spins and the experimental magnetic susceptibility, which 
allowed to determination of the EF directly from the experimental data. 

This paper is organized as follows: in the section |IT] we give a brief description of the particular system studied 
here. The following sections, third and fourth, contain the experimental results at zero applied field and a theoretical 
study of the entanglement evolution as a function of temperature and applied field respectively. In the last section, 
some comments and conclusions are drawn. 
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FIG. 1: The structure of the Na2Cu5Si40i4 compound, being the smaller circles the Oxygen atoms and the larger ones the 
Copper atoms, from two different views, (a) Side view, with an arrow indicating the staggered staking of the chains, (b) Top 
view of the chain. 



II. THE SPIN CHAIN DESCRIPTION 



The Copper atoms in the compound Na2Cu5Si40i4 have S — 1/2 and are separated in two groups, containing two 
and three atoms each, named dimer and trimer respectively. The whole structure of this compound is comprised of 
zig-zag clusters of copper and oxygen atoms then forming dimer-trimer sets of spins, as shown in Figure ([1]). There is 
an indirect exchange interaction between the spins, through the electronic oxygen clouds, and is of short range order, 
i.e., mainly between the first neighbors. As previously shown [3^ . the three spins that form the trimer are coupled 
antiferromagnctically with each other, whereas the two atoms of the dimer are coupled ferromagnetically. In addition, 
the two sets of spins, dimer and trimer, interact antiferromagnctically with each other. Therefore, labeling the spins 
according to the Figure ([2|), the magnetism of the this system can be described by the Hamiltonian: 



H = -MSi ■ 5*2 + 52 • 53) - J2(5; • Sb) - JsiSi ■ S5) - gfiBH ■ S (1) 

where g is the Lande factor, /x^ is the Bohr magneton, H is an external magnetic field, Sa = 84 + 85 is the total 
spin of the dimer, 8b = 51 + 52 + 53 is the total spin of the trimer and 8 = Sa + Sb stands for the total spin of 
the cluster. The values of the exchange integrals were determined experimentally (3^ . and they were found to be 
Ji = —224.9 K, J3 = 40.22 K and J2 = —8.01 K, for intra-trimer, intra-dimer and dimer-trimer, respectively. 

At finite temperatures T, the thermal equilibrium quantum state of the system is described by the density matrix 
p{T) = cxjp{—Ti./kBT)/Z where Z is the partition function. From p{T) it is possible to calculate thcrmodynamical 
quantities, such as magnetic susceptibility. Since the Hamiltonian Ti. commutes with the spin component along the 
z-direction 5^, one can show that the magnetic susceptibility along a given direction a can be written as [33 |: 



X"(T) = (^Y. (5f 5;) - I^Y. j (2) 

On Figure ([3]), the theoretical prediction of Equation Q is compared to the experimental magnetic susceptibility 
measured using a conventional superconducting quantum interference device (SQUID) magnetometer with an applied 
field of 100 Oe as a function of temperature. As can be seen from the figure, there is a monotonic increase of 
susceptibility down to 8 K, where there is a sharp drop associated to a transition to a 3D ground state [32|. The 
good agreement between the calculations and the experimental data above 8 K shows the faithful of the dimer-trimer 
cluster model described by Hamiltonian ([T]). 
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FIG. 2: (Color online) Schematic representation of the dimer-trimer cluster, and the respective interactions (J^) between the 
Cu atoms, according to Equation ([l}. 
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FIG. 3: Magnetic susceptibility as a function of temperature with an applied field of 100 Oe. The points (•) are the experimental 
results and the solid line is the theoretical prediction, based on Equation ([2}. 



III. ENTANGLEMENT FOR H = 



In this section, two quantities are used to investigate the presence of entanglement: an Entanglement Witness and 
the Entanglement of Formation, both obtained from the magnetic susceptibility. This study was carried out for the 
whole system and also for subsystems comprised by two and three spins. 



A. Magnetic Susceptibility Witness 



The concept of Entanglement Witness was first introduced by Horodecki et al. (31| , where it was defined that an 
EW is an observable which is capable to identify if a system is in an entangled state. However, in most of cases, an 
EW can not quantify the amount of entanglement present in the system and it gives only a sufficient condition for 
the presence of entanglement. In other words, an EW can only tell if, for some condition, a given system is in an 
entangled state, but if this condition is violated, we can not state with absolute certainty that the system's state is 
separable. 

Up to now, various thermodynamical quantities, such as magnetization [l^, internal energy p^ . [l^ . [Tsl. [TgI. [itI . [l9l . [20j 
and heat capacity [1^ , have been proposed as an EW, but these depend on the complete knowledge of the specific 
model describing the system and usually are not directly measurable quantities. However, recently, it has been 
demonstrated, by Wiesniak et al. [s^], that the magnetic susceptibility can be used as an EW for a large class of 
systems. This witness can be measured directly and can be applicable, in principle, without the full knowledge of the 
model Hamiltonian. If a system which has the symmetry [H, S^] = is in an entangled state, then the average of the 
magnetic susceptibility x'^^^'iT) measured along the three orthogonal axis satisfies the relation [s^ : 
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FIG. 4: Entanglement Witness for the total system EW{5) circles (•), for the trimer EW{3) triangles (a) and for the dimer 
EW{2) stars (*). 

Where N is the number of spin— particles. From Equation ([3]), we can define the Entanglement Witness as being: 



EW{N) 



(4) 



Thus, systems with EW{N) < are in an entangled state. However, it must be emphasized that the condition 
EW{N) > does not necessarily implies in separability. With EW{5) obtained from the measured magnetic suscep- 
tibility, we can determine the presence of entanglement in our 5 spins system. However it is also interesting to study 
if the subsystems are in an entangled state. In order to analyze the entanglement between spins in a subsystem, it is 
necessary to derive the contribution of the susceptibility Xsub(T) due to only the spins that matter. Theoretically, this 
can be done according to Equation but using the reduced density matrix Psub{T), which represents the density 
matrix of the subsystem. The Psub{T) can be obtained by using the partial trace operation 0], which sums over all 
the possible states of the spins in the system, except those belonging to the subsystem. From the numerical calcula- 
tion of the total susceptibility and subsystem's contribution, we can define the ratio Rl^TiJ") = X*t&°(^)/x*'"^°(^) 
which represents the fraction of the subsystem's contribution to the total magnetic susceptibility. With this quantity 
at hand, which can only be calculated with the knowledge of the Hamiltonian, it is possible to extract separately, 
from the experimental data, that contains the contribution of all spins, the magnetic susceptibility of the trimer 
XTrm = Rfr\"{T) X r^^T), dimer XdL^T) = R%Z{T) x r^^T) and so on. 

In Figure we show the experimental data, as can be seen, the EW{^) extracted from the total magnetic 
susceptibility is negative for any temperature below ~ 110 K, showing the presence of entanglement in the system. 
Furthermore, the fact that EW{?>) for the trimer is negative below ^ 240 K and EW{2) for the dimer is always 
positive suggests that entanglement only occurs between spins inside the trimer. 



B. Entanglement of Formation 



The Entanglement of Formation |29l . l30l | , which can be seen as the amount of quantum resources needed to create 
a given entangled state, is one of the most used for entanglement measurement. This was proposed to quantify the 
entanglement of a bipartite system. Unfortunately, as many others proposals, the EE is extremely difficult to calculate 
in general. However, in the special case of two spin 1/2 particles, recently an analytical expression was derived by 
Wootters who showed that the EF of any density matrix {p) of two spin 1/2 particles is given by [30| : 



EE ^ ~xlog2{x) - {I- x)log2{l- x) (5) 

where x = {l + yr^^C2)/2, being C the concurrence, defined as max(0, — -v/A^ — ~ ^f^) and A's are the 
eigenvalues of i? = pay ® ayp*ay ® Uy, labeled in decreasing order. The degree of entanglement varies from to 1, a 
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pair of spins is considered to be in a maximally entangled state if EF = 1 and separable when EF is equal to zero, 
for any other values the state of the spins is said to be partially entangled. 

Since the Hamiltonian ([!]) commutes with 5^, we can write the reduced density matrix of the spins located in the 
sites i and j as [H, HI] : 



f u+ 

uji z* 

z 

V u- 



being u± = (1 ± 2(5f + 5|) + 4(5f S'|))/4 and z = {S, 
density matrix can be written as: 



{S 



(6) 



- i{SfS^) - i{SfSJ). The concurrence of such 



C = 2max(0, \z\ - Vu+u~) (7) 

Now, exploring the fact that the system is isotropic at zero magnetic field, it is possible to write (SfSj) ~ {Sf S!-) = 
(S^Sj) = Gij/i and (SfS^) — (SfS^). Then, one can rewrite the concurrence between spins i and j in terms of 
the correlations functions, such as Cij{T) = | max(0, 2|Gy | — Gy — |). It is straightforward to verify from ^ that 
Xi^^{T) = 2(g/iB)^(l/4 + Gij/3)/kBT and thus the concurrence becomes: 



-xiriT?j (8) 

The Equation ([8]) relates the concurrence of the pair i ~ j, hence the EF, to the magnetic susceptibility Xij^{T), 
which can be obtained from the experimental data as explained before. It is interesting to note that the Equation ([5]) 
leads to the concurrence found by Asoudeh and Karimipour @ for mean-field clusters. 

In Figure ^ we show the experimental EF obtained from the measured magnetic susceptibility and the theoretical 
EF calculated with the reduced density matrix Pij{T) for the pairs inside the trimer (1 — 2,2 — 3 and 1 — 3). We 
can see that there is entanglement only between the pairs of spins 1 — 2 and 2 — 3 which persists up to a critical 
temperature ~ 200 K. The Entanglement of Formation for other pairs are always null and they are not shown. 
These results confirm that only the spins in the trimer are entangled, being those in the dimer in a separable state. A 
interesting feature is that the Entanglement Witness for the trimer EW{?>) gives ^ 240 while the obtained from 
the Entanglement of Formation is ~ 200 K. Since the EF can only see bipartite entanglement, this result indicates 
that pairwise entanglement is not the sole entanglement in the system and a tripartite entanglement is also present 
in the trimer. 

It is also interesting to ask if this multipartite entanglement is genuine or not. Genuine multipartite entanglement 
can be identified by the criterion described in (isj . Applied to the trimer system, the criterion states that if the 
expression {Utt;) < Ji{l + \/5)/4 holds, where {HTri) is the mean value of the energy corresponding just to the 
trimer part, then the tripartite entanglement is indeed genuine. We have numerically calculated, from the reduced 
density matrix of the trimer, the threshold temperature below which genuine tripartite entanglement exists. We 
found that the threshold temperature is ~ 108 K. This result suggests that the compound has a genuine tripartite 
entanglement for a finite temperature. Therefore, this genuine entanglement docs not extend up to 240 K. It is 
also important to emphasize the fact that entanglement is confined in the trimer and thus, it is not a macroscopic 
entanglement, but rather a thermal entanglement. 
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IV. THE EFFECT OF THE MAGNETIC FIELD 



In order to investigate the effect of the application of a magnetic field, the EF was calculated for each pair of 
spins. In Figures ^ and ([7]) we show the EF for the 3 pairs of spins in the trimer. The entanglement of the pair 
1 — 2 behaves as the same way as of the pair 2 — 3, the application of a low magnetic field increases the amount of 
entanglement of both pairs at low temperatures and also creates a small entanglement between the pair 1 — 3, which 
was zero at zero field. Similar behavior was also encountered in other theoretical works 0, and can be understood 
in terms of changes in the ground state. At iJ = 0, the ground state of the system is degenerated, being a statistical 
mixture of two states. This mixture has only entanglement between pairs 1 — 2 and 2 — 3. The application of an 
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FIG. 5: The experimentally determined Entanglement of Formation, the circles, for the pairs inside the trimer 1 — 2 (a), 2 — 3 
(b) and 1 — 3 (c). The solid lines are the theoretical prediction. The small discrepancy between theory and experiment at low 
temperatures is associated to a 3D transition [3^ . 
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FIG. 6: (Left) Entanglement of Formation (EF) for the pairs 1 — 2 and 2 — 3, as a function of temperature and applied field in 
the H-T plane. (Right) EF for selected temperatures, as a function of applied magnetic field . 



external magnetic field changes the energy eigenvalues leading to a different ground state which has a different degree 
of entanglement. 

For a high enough magnetic field, the ground sate becomes the state |tTTTT)i which correspond to all spins aligned 
parallel with the external field. At low temperatures, where the ground state is very populated, EF vanishes suddenly 
near ^ 2200 Oe due to the sudden change in the ground state to ITTTTT)- However for higher temperatures, there are 
many states populated and the decreasing of entanglement is slow. 
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FIG. 7: (Left) Entanglement of Formation (EF) for the pair 1 — 3, as a function of temperature and applied field in the H-T 
plane. (Right) EF for selected temperatures, as a function of applied magnetic field. 

Prom the H — T diagram showed in Figures ^ and (O, we can also observe that the entanglement of pairs do 
not occurs above the critical temperature ^ 200 K and above the critical field H'^ ~ 3000 Oe. A interesting 
feature is that the critical field H'^ is much higher than the available field intensities that usually can be create in 
laboratories, which indicates that the entanglement in this system can not be destroyed by common magnetic field for 
temperatures below 200 K, providing that high magnetic fields do not induce any structural changes in the compound. 
The Entanglement of Formation for others pairs are always zero at any point in the H — T plane and they are not 
shown. 



V. CONCLUSIONS 

In summary, we have successfully established the presence of thermal entanglement in the compound Na2Cu5Si40i4. 
From the obtained results, we could conclude that entanglement is strong in this system and disappears only at high 
temperatures and high magnetic fields. At zero magnetic field we found experimentally that pairwise entanglement 
exists below ~ 200 K and there is tripartite entanglement below ^ 240 K, which shows that multipartite entanglement 
is stronger than bipartite entanglement in our case. Same feature was also verified experimentally in the compound 
Na2V3 07 [13]. The entanglement observed here is confined into small cluster of spin and thus it is not a macroscopic 
entanglement. We have also presented a theoretical study of the entanglement evolution as a function of applied field 
and temperature, showing that magnetic field can increase the degree of entanglement in the pairs 1 — 2, 2 — 3 and 
1 - 3. 

We emphasize that some results of this paper is based on the validity of the dimer-trimer model, however this model 
is shown to be a good model for the present system (see Figure (|3])). Furthermore, since the Entanglement Witness 
EW{b) does not require any assumption about the model, or the explicit values of the exchange-coupling parameters, 
if the model is not completely correct, the main conclusion does not change, i.e the presence of thermal entanglement 
in the compound Na2Cu5Si40i4. 
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